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Abstract 

The computation of the one-loop effective action in a radially symmetric background can be 
reduced to a sum over partial- wave contributions, each of which is the logarithm of an appro- 



. priate one-dimensional radial determinant. While these individual radial determinants can be 

o , 

\G evaluated simply and efficiently using the Gel'fand-Yaglom method, the sum over all partial-wave 

o 

' contributions diverges. A renormalization procedure is needed to unambiguously define the finite 

"t^ '. 

renormalized effective action. Here we use a combination of the Schwinger proper-time method, 

<D '. 

r~| , and a resummed uniform DeWitt expansion. This provides a more elegant technique for extracting 

the large partial-wave contribution, compared to the higher order radial WKB approach which 
^ • had been used in previous work. We illustrate the general method with a complete analysis of the 

scalar one-loop effective action in a class of radially separable SU(2) Yang-Mills background fields. 
We also show that this method can be applied to the case where the background gauge fields have 
asymptotic limits appropriate to uniform field strengths, such as for example in the Minkowski 
solution, which describes an instanton immersed in a constant background. Detailed numerical 
results will be presented in a sequel. 
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I. INTRODUCTION 



In the study of quantum field theories, one is often led to consider the one-loop effective 
action in nontrivial background fields. While the renormalization counterterm structure of 
one-loop effective actions can readily be exhibited for general backgrounds, the explicit eval- 
uation of its full finite part in an interesting specific background still constitutes a highly 
nontrivial problem. For gauge theories, explicit analytic results are known only for very 
special backgrounds: the Euler-Heisenberg effective action for a background with constant 



strength 



, its generalization to a nonabelian covariantly constant back- 



For applications 



abelian fie 

ground J, |5|, |6|, l7| , and a special solvable abelian background 
in both continuum and lattice field theory, one would like to enlarge this set of backgrounds 
for which we have accurate computations of the finite renormalized effective action. 

In a series of recent publications (together with Hyunsoo Min) 2|; we presented a new 
method for computing the renormalized one-loop effective action in a radially symmetric 
nonabelian background, and used it to evaluate explicitly the QCD single-instanton deter- 
minant for arbitrary quark mass values. The related computation in the massless limit was 
performed in a classic paper of 't Hooft while the heavy quark mass limit was studied 



in Q 



15l |. The new method in [12] works for any quark mass, not relying on small or 



large mass expansions, and the result interpolates smoothly and precisely between these 
two extremes. In this paper we present the general formalism, and we introduce a simpli- 
fied analysis based on a uniform Schwinger-DeWitt expansion, replacing the higher order 



radial WKB analysis used in 



12 



llGl l. This approach has also been used to evaluate the 

n 

exact determinant prefactor in false vacuum decay (or nucleation) [17|. We also note that 
a related method has recently been applied to the two-dimensional chiral Higgs model [3| ■ 
Finally, an alternative derivation, using zeta functions rather than a partial-wave cutoff, of 
the determinant of a radially symmetric Schrodinger operator has been given in [igj ]. 

The starting idea is very simple. If the background field is radially symmetric, the effective 
action F can be expressed formally in terms of one-dimensional functional determinants of 
radial differential operators for various partial waves. Explicitly, writing J for all quantum 
numbers specific to a given partial wave, it has the general structure 

detij-ij + w?) 



CO 



~ ^In 

j=o 
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Here Tij, the radial Schrodinger-type differential operator for the J-th partial wave, contains 
a nontrivial (background-dependent) potential term, while Ti^j'^^ is the corresponding free 
operator. In general there will also be appropriate degeneracy factors in the sum. For a 
given partial wave J, the individual determinant is finite once we divide by the corresponding 
'free' contribution. These finite one- dimensional determinants can be evaluated easily using 
the Gel'fand-Yaglom method which reduces the computation to 

a trivial (numerical) integration with initial value boundary conditions. 

The non-trivial aspect of this approach is that the sum over partial waves J diverges (in 
spacetime dimension d > 2). This is of course because the formal sum in p. HI neglects 
renormalization. We solve this problem by introducing a partial- wave cutoff Jl, and we 
isolate the divergence of the sum in a form that can be absorbed via renormalization. Specific 
details are presented in the body of this paper, but the general structure is that we write 
the sum as 

^ fdet{Hj + m^)\ / det(7^J + m2) \ ^ ^ .^ 

r = y In ^ ,/ ^ + y In ^ ,/ ^ + counterterm . 1.2 

The finite, renormalized effective action is then evaluated as follows: 

• The first sum is evaluated numerically using the Gel-fand Yaglom method. 

• The second sum is evaluated analytically in the large Jl limit, after regulating the 
determinants. This step uses our uniform Schwinger-DeWitt expansion and Euler- 
Maclaurin summation. 

• The analytic large Jl behavior of the regulated determinants leads to the correct 
renormalization counterterm, and moreover cancels exactly the numerical divergences 
of the first sum as — > oo. This produces a finite renormalized answer. 

The technically difficult part of the computation is the analytic computation of the large 
Jl behavior of the second sum in (11.21). This was achieved in using second-order radial 
WKB, based on Dunham's formula 27|. While this is very general, it can be quite cum- 
bersome for complicated background fields. In this paper we present a simpler method to 
implement this part of the computation. The analysis reduces to simple algebraic manipu- 
lations, and can be more readily generalized. With this new approach we can now evaluate 
exactly the finite renormalized effective action for a very general class of radially-symmetric 



backgrounds. In fact, with a background field involving an unspecified radial function, the 
large partial-wave contribution to the renormalized effective action can now be evaluated 
explicitly. This will be important to discuss the background-field-dependence of the effective 
action and also to test various approximation schemes. 

This paper is organized as follows. In Section IH] we define the renormalized one-loop 
effective action for the case of a scalar field in a class of spherically-symmetric Yang-Mills 
background fields, assuming four dimensional Euclidean spacetime. Its Schwinger proper- 
time representation [2] is given. The full amplitude is then expressed using partial wave 
amplitudes, and we also elaborate here on the role of the two kinds of proper-time Green 
functions, one related to the quadratic differential operator given in 4-D spacetime and 
the other for the radial quadratic differential operator. In Section lllll we explain how the 
large-/ partial wave contribution to the full effective action can be evaluated explicitly using 
a generalized DeWitt WKB expansion for the radial proper-time Green function (i.e., the 
y-expansion) and the Euler-Maclaurin summation method. This allows us to present the 
renormalized one-loop effective action (in the class of spherically symmetric backgrounds) in 
a form amenable to direct numerical analysis. In Section HVl we show that our y-expansion 
formula can be applied to the calculation of the large partial-wave contribution even when 
background gauge fields do not fall off at large distance but approach those of uniform 
field strength. Also given here is the exact partial-wave-based treatment of the effective 
action in the background corresponding to strictly uniform self-dual field strengths. In 
Section |3 we conclude with some relevant discussions and comments. There are several 
appendices which contain supplementary materials and some technical details. In Appendix 
^the explicit form of the free radial proper-time Green function in n spacetime dimension 
is considered. In Appendix El we study the coefficient functions in the y-expansion when 
the potential is matrix- valued. Appendix contains a brief account of the Euler-Maclaurin 
summation formula, and some explicit results obtained using this formula in connection with 
our problem. In a sequel we will address matrix-valued problems in more detail, and present 
detailed numerical results for the general radial cases for which the formalism is developed 
in this current paper. 
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II. EFFECTIVE ACTION IN RADIALLY SYMMETRIC BACKGROUNDS 



Renormalized Effective Action 



We choose for our field theory model an SU(2) Euclidean Yang-Mills theory with a com- 
plex scalar matter field (in the fundamental representation), in four dimensional spacetime. 
As far as our general methodology is concerned, the model choice is not crucial; but, by 
choosing this case, we are able to crosscheck readily the findings of the present work against 
those of Refs. llG^ . The case with Dirac fields is quite similar if one works with the 
squared Dirac operator. Also, by considering a gauge theory (rather than the much simpler 
scalar field theory), we can demonstrate the gauge invariance of our calculational method 
for the renormalized effective action. 

Consider a generic Yang-Mills background : A^{x) = 74^(x)^ (/x = 1,2,3,4, a = 1,2,3; 
r's denote 2x2 Pauli matrices). The Pauli-Villars regularized one- loop effective action 
associated with scalar field fluctuations can be represented by 

'det{-D^ + m^) det(-92 + A^) 



V tS^A; m) = In 



(2.1) 



_det(-92 + m2) det(-D2 + A2)_ 

where m is the scalar mass, A a heavy regulator mass, and the covariant Laplacian 
operator 

D^ = D^D^, , {D^ = d^-tA^{^)). (2.2) 
The Schwinger proper-time representation for the form ()2.H) is 



°° ds 
s 



TA{A;m) = 

with F{s) given by {s is the proper-time variable) 



F{s) 



(2.3) 



F{s) = J rf^x tr (x|{e 
The proper-time Green's function 

A(x, x'; s) = (x|e^ 



(2.4) 



(2.5) 



admits an asymptotic expansion, the DeWitt (or heat kernel) expansion 28|, |29| : 

1 



(47rs)^ 



,-|x-x'|''/4s 



s"a„(x, x' 

.n=0 



for s — s> -|- . 



(2.6) 



The expansion coefficients, a„(x, x') (n = 0, 1,2, ■ ■ ■), and especially the coincidence limits 
a„(x, x) of the first few terms, can be found most simply using recurrence relations satisfied 
by the a„(x, x')'s. The divergence structure of T\{A; m ) as A oo is governed by the values 
of tr ai(x, x) and tr 02 (x, x), and in our case we have j2^ 



trai(x, x) = 0, tra2(x, x) 



--tr[F^,(x)F^,(x)] 



(2.7) 



where F^^, = F^y\ = i[Dfj_,Dy\ is the field strength. Then the renormalized one-loop 
effective action in the minimal subtraction scheme is defined as 



lim 

A— >oo 



V^{A-m) - ^^In J rf^xtr[F^,(x)F^,(x)] 



(2.8) 



where n is the renormalization scale. 



B. Radial Backgrounds 

It is a very difficult problem to explicitly evaluate this renormalized effective action ()2.8p . 
For a generic background gauge field y4^(x), there is currently no known method leading 
to an exact evaluation of the one-loop effective action. On the other hand, there are many 
interesting physical applications {e.g., vortices, monopoles, instantons, . . . ) where the gauge 
background is radially symmetric. In this paper we show that this radial symmetry is strong 
enough to permit the computation of the renormalized effective action ()2.8p . 

A large class of such radial backgrounds is covered by the ansatz form: 

A^(x) = 2r]^^aXuf{r)— + 2ri^^:iXyg{r)— , (r = |x| = ^x^x^) (2.9) 

where the radial functions /(r) and g{r) are left unspecified, and 77^,^0 (a = 1,2,3) denote 
the standard 't Hooft symbols [13]. With ^^^(x) of the form ()2.9|) . the covariant Laplacian 
operator —D^ becomes (here ^ = T^) 

- -D^ = -d^d^ + 4:if{r)r]^i,aTaXudf, + Aig{r)7]^^^T^Xud^ + 4:f{rY7]^yail^i\bTaTbXuXx 

+4:g{rfri^^^ri^xzTixyXx + ^f{r)g{r)ri^yarit^xzTaTzXyXx • (2.10) 

We may then define [3| the operators La = —i^VpLvaXptdy (satisfying angular-momentum 
commutation relations [La, L}] = ieabcLc) and use the relations 

3 

Vfiua^tiXb = ^ab^uX + ^abcVuXc , TaTa = -1 , 

-^A = -|^-'|: + ^^^ , {L'^LaLa) (2.11) 
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to recast the expression ()2.1U|) as 

92 3d 



+r'{3f{rY + g{ry + 2f{r)gir)} . 



(2.12) 



Based on this form, we may associate an infinite number of partial-wave radial differential 
operators with the given system. We distinguish between three important cases. 

1. Case 1 : g{r) = 0, hut f{r) ^ 



Suppose that q{r) = 0, but /(r) 7^ 0. This is the form relevant to the instanton com- 



putation in 



13l |. Then A^(x) is given by the first piece only on the right hand side 



of ()2.9|) . Then, noting that there exists another set of angular-momentum-like operators 
fj^yaX^dy (satisfying [La, Lb] = and LaLa = LaLa = ) jlSj, partial waves can 



"2 



be specified by the quantum numbers Ji = {l,j,j3,h), where 
(L7 = /(/ + 1), / = 0,^,1,^,--- ; 
{•Py = j{j + 1), (with J, = L, + TJ, J 

{J^y = 33 = -j, -j + 1,- - ■,j ; 
(L,y = i, = -l,-l + i,...,L 



(2.13) 



The radial differential operator, representing —D^ in the given partial wave sector, thus 
assumes the form 

3- 



n 



Ji 



-D 



5(0+4/(r) 



j(j + l)-/(/ + l) 



+ 3rV(r) 



where 9^^;^ is the partial-wave form of the free Laplacian 9^9^ : 



(2.14) 



(2.15) 



2. Case 2: f{r) = 0, but g{r) / 

The system with /(r) = 0, but g{r) 7^ 0, is simpler: here, partial waves are specified by the 
quantum numbers J2 = (/, /s, ts, Is), where {L^y = 1^ = —I, — / + 1, ■ ■ ■ , / and (T3)' = = ±|. 
In this case, the radial differential operator becomes 



Hj, = -D 



-dl) + 9,g{r)ht^ + r^g{r) 



(2.16) 
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3. Case 3: both f{r) and g{r) nonvanishing 



The situation is somewhat more comphcated if both /(r) and g{r) are nonvanishing, since 
—D^ will then be nondiagonal in either basis considered above. This case can be treated by 
allowing partial wave sectors themselves to be finite-dimensional vector spaces. Explicitly, 
taking a partial wave specified by the quantum numbers J3 = (/, js, I3), we can represent the 
operator T3L3 according to 



/ 

21+1 4 21+1 



13 (1 I l^2 _ a2 J'i _ 1 2 2 

V 2«+l V " 2/ 2i+l 4 



T3L3 ^ h, if j3 = ±(/ + ^), (2.17) 

where the 2x2 matrix, appearing when |j3| 7^ / + |, is defined relative to the basis vectors 
|j = / ± |). This allows us to represent —D^ in the given partial wave by the (matrix) radial 
differential operator 

^ ^_j,2 ^ f + W^(0 + %.3) ,ifj3 = -/ + |,---,/-| 

^''''^ 1 -dfn + Wir) + Alf{r) + 4/^?(r) , if j, = ±(/ + |), 



where 



H^(r) = {3/(r)2 + g^r^ + 2/(r)(7(r)} , (2.19) 
^(ija) = — , (2.20) 



4//(r) + 8(^ - \)g{r) -^^^{l + \y - j| ^(r) 

-Hh^il + \?~3l9{r) -4(/ + l)/(r)-8(^ + i)^?(r) 



The Gel'fand-Yaglom method has a straightforward generalization [25[ to matrix- valued op- 
erators, so the numerical part of the computation follows as before. Such matrix-valued 
radial operators have in fact been considered in and also occur naturally when consid- 
ering fiuctuations of a Dirac-spinor matter field in a radially symmetric background. 

An interesting subclass of these radial backgrounds consists of those that are self-dual 
(or anti- self -dual). Such gauge fields satisfy automatically the classical Yang- Mills field 
equations, and as such they are of particular importance. With our potential form in (j2.9|) . 
such self-dual or anti-self-dual configurations are obtained if certain special functional forms 
are chosen for /(r) and g{r). Explicitly, for self-dual configurations, the following choices 
can be made: 

(i) /(r) = :pi^, g{r) = (i.e., = 'q^yaipr^T'') for a single instanton solution in the 
regular gauge; 



(ii) /(r) = 0, g{r) = = const, (i.e., = —rf^yaX^^'T"') for a uniform self-dual field 
strength background; 

(iii) /(r) = sinh[&(r2^p2)] ) gi''^) = & tanh[|(r^ + p^)] for the so-called Minkowski solution [s^ 
which describes a single instanton immersed in a uniform background. 

Note that the Minkowski solution reduces to the case (i) or (ii) in appropriate limits. Anti- 

2 

self-dual solutions may be obtained for the choice /(r) = ^2(r2_|_p2) ; and g{r) = (i.e., 

2 

^fj. = ^/Jt/a ^2(p2^p2) ); this corresponds to a single anti- instanton in the singular gauge. With 
any of these classical solutions chosen as the background, our discussion above tells us that 
the operator —D^ can be written in the partial-wave expanded form. In particular, in the 
case of the Minkowski solution for which both /(r) and g{r) are nonvanishing, the related 
partial-wave differential operator will take a 2 x 2 matrix form. 

But, in the following analysis, it will be sufficient to assume that our background po- 
tentials are just of the radial form ()2.9p — i.e., they do not have to satisfy classical field 
equations. 

C. Partial Wave Decomposition of Effective Action 

Taking advantage of this radial symmetry, we can make a partial wave decomposition in 

F{s) = Y.Fjis) , (2.21) 
J 

where 

Fj{s)= J^^drtT{Aj{r,r;s)-A'r{r,r;s)} . (2.22) 
The proper-time radial Green's function for each partial wave J is defined as 

Aj(r,r';s) = (r|e-"^'^ |r') , (2.23) 

in terms of the radial operator 

nj^^,njr'/' = -^,+VAr) . (2.24) 

Note that we have extracted a measure factor r^^"^ in writing H = r^^/'^Tir^/'^ . The form of 
the (possibly matrix valued) radial potential Vj(r) depends on the specific form of the gauge 



9 



field entering the covariant Laplacian operator. In each case, the effective radial partial-wave 
potential Vj(r) contains a centrifugal term, having the structure 

VAr) = ^^^^±^ + UAr) , (2.25) 

where / is the half-integer valued quantum number in ()2.13|) . 

The partial-wave-based representation of the regularized effective action is : 

TK{A-m) = -Y^j^ -(e-^'^-e-^'') dr tr {A,(r, r; s) - AS^-(r, r; s)} . (2.26) 

Then, for the explicit evaluation of the renormalized effective action given by ()2.8p . it is 



convenient to separate the partial wave sum into two parts [12|, 1^ : (i) the sum over partial 
waves with J < Jl (here, Jl is chosen such that it may refer to some large Z- value, I = L); 
and (ii) the remaining infinite sum involving all J > Jl terms. In the first contribution 
involving the finite J-sum, the regulator plays no role in the limit A — oo, and so may be 
removed from this sum. Based on this procedure, we can now write 

r,en(A;m) = rj<j^(A;m)+rj>jJA;m) , (2.27) 



with 



JZjJ'> s Jo ^ J 

Elndet(l^) , (2.28) 



J<Jl 

r°° ds 



rj,j,{A-m) = - E r ^ (e-^^ - e-^'^) Fj{s) 

In the right hand side of (IT^ . lndet( ^;^+_^J may well be replaced by lndet( ^4+_^2 ), 
the two being the same. In our model cases, we notice that not all quantum numbers in J are 
relevant for our radial Green function Aj(r, r'; s). In view of this, for the three possible forms 
of a. J considered above [see ()2.14|) . ()2.16|) and ()2.18p ]. we may express the decomposition 
formula ()2.21|) in the more explicit forms: 

Case 1 : F{s) = ^(2j + 1)(2/ + l)F(,,,)(s) 

10 



E (2/ + 1)(2/ + 2)[F(,,^.^,+ i)(.) + , (2.30a) 

z=o,^,i,--- 



Ca5eS:F(s)= ^ ^ ^ (2^ + WMa.ta)!^) 

1=0,^,1,- h=-lt:i=±^ 

Case 3: F{s) = ^ ^ (2/ + l)F(,,,3)(s) . 

^=o,|,i,-i3=-a+|) 



(2.30b) 



(2.30c) 



Note that the degeneracy factors here contain a common factor (2/ + 1) from the /3-sum. 
Here, as the notations of (j2.3(Ja|) - (j2.3(Jcp are used, the designation J < 01 J > may be 
identified with the appropriate division in the values of the quantum number /, i.e., J < Jl 
when I < L {L is some, arbitrarily chosen, large value) and J > Jl when / > L. For 
example, the low partial wave sum in ()2.28|) is expressed explicitly as: 



Case 1 : 



det (-1^5,^+1/2) +m2) det(-D2 + 



+ ln 



det(-95+,/2) + "^' 



Case 2 : 
Tj<j^{A,m) 

Case 3 : 



J-^ ^ det (-£>?,,, ^ 

^ ^ E (2^ + 1) In ■ ^ • 

Z=0 h=-lt3=±l/2 



det(— (9^;^ + m^) 



j<jj 



(A, m) 



L '+5 



det(-^L3)+"^') 



y y (2/ + i)in , 

i^..^^^) det(-55)+m2) 



(2.31a) 
(2.31b) 

(2.31c) 



The low partial wave contribution, rj<j^, in p.28|l . may be determined numerically. On the 
other hand, the large partial wave contribution, rj>j^, in (|2.29p . is calculated analytically 
for large L, to the desired accuracy in powers of 1/L. 



D. Low Partial Wave Contribution 



To evaluate the lowpartia 



nique 



21 



22 



23 



25 



wave contribution F 



J<Jl 



we use the Gel'fand-Yaglom tech- 



which can be summarized as follows. Suppose M.i and 
M.2 denote two second-order radial differential operators on the interval r G [0, 00). Then 
the ratio of the determinants is given by 



det R-^oo\^2{R)^ 



(2.32) 
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where $i(r) {i = 1,2) satisfy the initial value problems : 

Mi^i{r) = ; <l>i(r)~r2' as r^O . (2.33) 

Here / is the index in the centrifugal term in ()2.25|) . Since an initial value problem is ex- 
tremely simple to solve numerically, this provides an efficient calculational method for the 
individual radial determinants. Here we take M.i = Hj + m^, and = 'Hj'^'^ + m?. 
Thus, ^2{j') is in fact known analytically. Then better numerical results are obtained 
by considering directly the initial value problem with the second-order differential equation 
derived for the ratio function S{r) = ln($i(r)/$2('"))- This method has been implemented 
successfully in for the instanton determinant computation, and for the false vacuum 



decay problem in both flat 



17| and curved 



spacetime. Furthermore, this method of 



calculating radial determinants can be generalized to the case when the second-order dif- 

I I 

ferential operator Ai in question contains a matrix-type potential 25], as in our Case 3. 
Explicit numerical results for the three radial cases discussed above will be presented in the 
sequel. 

III. LARGE PARTIAL- WAVE CONTRIBUTIONS AND RENORMALIZATION 

The large partial wave contribution cannot be evaluated numerically because of the need 
to remove the heavy mass regulator, A. Instead we compute analytically the large L behavior 
of Tj^j^. To determine the large partial- wave contribution Tj^j^ (given by (|2.29|)), o ne needs 



the large-/ (i.e., I ^ L) behavior of the function Fj{s). To that end, in Refs. we used 

the scattering phase shift representation of Fj{s), and then the radial WKB approximation 
up to second order. Here we present a much simpler approach, introducing a new "uniform" 
DeWitt expansion for the radial proper-time Green's function Aj(r, r;s), which remains 
valid when / becomes large. As we shall see, this new approach gives rise to results in 
complete agreement with those from our earlier method, with much less labor. 

A. Uniform DeWitt Expansion 

In the presence of the effective radial potential, V{r) = '^'^'''^2^'''^ + V(r), we seek a large-/ 
asymptotic representation of the related proper-time Green's function A(r, r'; s): 

{ds - dl + \/(r))A(r, r'; s) = , (for s > 0) (3.1a) 
12 



s ^ 0+ : A(r,r';s) — >S{r-r') . (3.1b) 

We take V(r) to be a typical smooth potential. Now, as / becomes very large, the presence 
of the large centrifugal potential, ^^^'^2"*^^ , has the consequence that our Green's function 
A(r, r'; s) acquires a totally negligible amplitude for s ^ A/P (where A is an 0(1) constant) . 
Because of this property, given a certain quantity which involves the integral of this function 
over s [such as, for example, rj>j^(A; m), given by (|2.29p ]. it will suffice to use an accurate 
representation of A(r, r';s) for s satisfying the condition < sl^ ^ 0{1), and having the 
property that it becomes exponentially small for s ^ A/l'^. Although only small-s values 
are relevant, the usual small-s DeWitt expansion [the one-dimensional analogue of ()2.6p ] 
cannot serve this purpose since it fails to account for the effect of the large centrifugal 
potential term. There is a conflict between the small s limit and the large / limit. To see 
this problem more clearly, consider the behavior of the function A(r, r'; s) with V(r) set to 
zero. For this free case, denoted A^'^'^'^^r, r'; s), we have a closed-form expression in general n 
spacetime dimensions (see Appendix^. As / becomes large, this function admits a uniform 
approximation of the form (see ()A10|) ) 



A'''-{r,r';s) ^ ^=e-^ {l + 0{s^)} , (3.2) 

valid as long as s is such that < sZ^ ^ 0{1). The naive small-s DeWitt expansion is 

4i(i+l)+| 

not adequate for our purpose since it effectively replaces the exponential factor e ^r' ^ 
(which can be 0(1) for s ~ ||) by the first few terms of its Taylor series in s. 

To obtain the desired large-/ expansion of our radial proper-time Green's function, it is 
convenient to set 

V{r) = l'^U{r) , (3.3) 
(so that U{r) remains finite as / 00), and introduce the rescaled proper-time variable 

t = fs . (3.4) 

Now the situation for the large-/ limit of A(r, r'; ^) is actually the same as that appropriate 
to the so-called -^-expansion of the proper-time Green function considered previously (for a 
different purpose) in Ref. 3^, identifying A with P. Thus, based on the result of |3^, we 
may immediately write the y-expansion of A(r, r'; ^), having the structure 

A(r,r';-) = ^e-^^<|E&.(^,^';i)(7^) \ , (3-5) 

<k=0 
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with suitable coefficient functions bk{r,r';t) which are regular near r = r' . Inserting this 
form in ()3.1a|) . we see that the coefficient functions bk{r,r';t) must satisfy: 

I'Uir) Y: h{r, r';t)(-) + ^ dMr, r'; t) 



fc=0 ^ fc=0 



Y.dMry;t)[-,) -T.9%{ry;t)[-^) =0. (3.6) 



^ k=0 ' k=0 



We can here regard f/(r) to be strictly of order {p)^, i.e., disregard the fact that it might 
contain terms with j-suppression, to simplify the presentation of our result. Then, 
gives rise to recurrence relations satisfied by the coefficient functions bk{r,r']t): 



0{f) : f/(r)6o(r,r';t) + 9t6o(r,r';t) + ^^9A(r,r';t) = 0, (3.7a) 



0{l'-'') ■ U{r)bk{r, r'; t) + dMr, r'; t) + —^dMr, r'; t) 

-d%.iir, r'; t) = 0, (A; = 1, 2, 3, ■ ■ ■) (3.7b) 

Further, because of the boundary condition ()3.1b|) . we must have bk{r,r';t = 0) = 5^0- 

To simplify the analysis of the recurrence relations, we may introduce a new variable u 
(instead of r) by setting r = r' + tu and define a new set of functions, 

h{u,r',t) =e^Ir''''^^^''^^'"bk{r' + tu,r';t). (3.8) 

[Here we have restricted our attention to the case when U{r) is not a matrix-valued potential. 
The case with a matrix- valued potential is discussed in Appendix Then the above 
recurrence relations can be recast as 

0{f) : ^bo{u,r';t) = 0, (3.9a) 
Oil'^'") : ^6,(«,r';t) = i||^6,_i(«,r';t)-2(?'(«,r';t)|-6fc_i(n,r';^ 



+ 



g'{u, r'- tf - g"{u, r'; t)] h-i{u, r'; t)} , (fc = 1, 2, ■ ■ ■) (3.9b) 



where g{u,r';t) = ^ J^!^^"" U{w)dw, g'{u,r';t) = ^g{u,r';t), and g"{u,r';t) = -^g{u,r';t). 
Since bQ{r,r';t = 0) = 1, and so bQ{u,r';t = 0) = 1, we now immediately conclude from 
fj3.9a|) that bo{u,r';t) = 1 for any t > 0. This in turn tells us that 

6o(r,r';t) = e-^ (3.10) 
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As a check, we note that if we choose the form U{r) = '^'^'^^2^"^ (appropriate to the free case 
with V(r) = 0), then ()3.10|) reduces to 60 (^"5 r']t) = e ^ f^, producing the correct expo- 
nential factor in ()3.2|) . related to the centrifugal potential term. Clearly, in the coincidence 
limit of r' = r, we have 

6o(r,r;t) = e-*^«. (3.11) 
Using bo{u, r'; t) = 1 in the k = 1 case of ()3.9b|l . we obtain 

^61 (u, r'- 1) = i [g'{u, T'-tf~ g"iu, r'; t) } . (3. 12) 

Then, to find the coincidence limit of 6i(r, r';t), i.e., the expression for r' = r, we may set 
M = in ()3.12|) (together with the easily obtained expressions g'{0,r';t) = |t^?7'(r) and 
g"{0,r';t) = lt^U"{r')) to obtain 

p,iu = 0,r;t) = \t'U'irr - V'(r). (3.13) 

This immediately leads to the expression 

Higher-order coefficients can be found similarly; for instance, for b2{r,r]t) we find 

b,{r,r;t) = e"*^^''^ " ^t'f/'W^f/"(r) 

+^tV(r)2 + ltV(r)f/(3)(r)-lt=^^W(r)} . (3.15) 

We can now exhibit the desired y-expansion structure for our radial proper-time Green's 
function in the coincidence limit. Returning to the notations using V{r) and s, it takes, 
based on the results of ()3.1ip and ()3.14p . the following form 

A(r, r; .) = ^e'^^W {l + {^s'V'{rf - ^.V"(r)) + O (^) } . (3.16) 

[See Appendix ^for the explicit verification that this gives rise to a correct large-/ expansion 
for A^'^'^^{r, r; s)]. In effect we have resummed all non-derivative terms in the standard DeWitt 
expansion. These non- derivative terms are all of the form {sV)^ for some k, and recalling 
that V depends quadratically on the partial wave index /, we see that all these terms are 
of 0(1) for s/^ ~ 0{1). On the other hand, the remaining terms in the expansion ()3.16|) go 



3 
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like s^l^ ~ s(s/^)^, and s'^P ~ ^(s/^), etc, and so remain small in this uniform limit of small 
s and large 1. 

This modified expansion may be used even when V{r) contains, apart from 0(/^) terms, 
some subleading terms as with the case V{r) = l^U{r) + lT{r) + Q(r); in this case, one 
can also generate an equally-valid y-expansion starting from the form (j3.16p by having the 
exponential of the subleading terms, i.e., e"'^^^^^''-*"'''^'''^^^ expanded (partly or wholly) in 
powers of j. This can be justified when r is restricted to the range in which T(r) and 
Q{r) remain bounded. This trivial rearrangement can, in fact, be incorporated within our 
y-expansion ansatz by allowing the power series development in the ansatz to have 

also odd-power terms in y. 



B. Explicit Large Partial Wave Contributions 

We may use the form with the formulas ()2.22|1 and ()2.29j) . to determine explicitly 

-wave contribution to the effective action. To facilitate this calculation, 



the large partia 

n , 

we follow Ref. [Ij, ll6| by trading the regulator mass A for a dimensional regularization 
parameter e. This is achieved by demanding that 

_ f-rn^s _ g-A^A p^^^ _ _ r ^e-™VF(s). (3.17) 

Jo s ^ ' Jo s 

Since F{s) = (finite constant) + 0{s) for small s, we then see that ()3.17p requires 

-ln(^^^ +0(^) =-^ + (7 + 21nm) + 0(6), (3.18) 
where 7 = 0.5772... is Euler's constant. Thus the relation between e and A is given by 

7 + m A^ 

Note that this is only to simplify our calculations; all the s-integrations appearing below 
can also be carried out within the original Pauli-Villars regularization framework. 

With this preparation, we now proceed to the calculation of rj>j^(y4; m) for our Case 1 
and Case 2. Case 3 will be considered in the sequel. 
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1. Case 1 



With Tijj given in ()2.14|) . we have the radial potential 



= ^^^^ + 4/(r) { j{j + 1) - /(/ + 1) - 



3rV(r) 



(3.20) 



which may be used in ()3.16p . Then, from ()2.22|) . Fj{s) for large enough / will follow. 
Further, if we represent F{s) by the form ()2.30a|l and use the correspondence ()3.19|) . it is 
possible to express the first part of ()2.29j) [the contribution to rj>j^(y4; m) other than the 
renormalization counterterm] as 



J>Jl 



{A;m) 



oo roo 

dr / ds 
Jo 



^-mH^e J2 (2/+l)(2/ + 2) 



{A(^^^.=;+i)(r,r; s) + A(;+i_ .^;)(r, r; s) - AJ^f (r,r; s) 



A;rA)(r,r;.)} 



(3.21) 



where we placed the s-integral before the r-integral, in order to give an explicit result for 
Tjyj^{A; m) for general fields. Using ()3.16p in ()3.21|) with V(i,j)(r) given by ()3.20p . the right 
hand side of 1)3.211) can be expressed as 



^ J>Jl 



[A] m) 



oo /-oo 

dr / ds 
Jq 



V4 



-sV, 



1 + 



g-m s^e J2 (2/ + l)(2/ + 2) 



+ 1^ 



sV, 



1 + 



12" 'C+i'O 



-S V(j,l n 



"(r) 



O 



/4 



-e-^^r(^) (^1 + 



12^ 



+ 



/4 



where V(;)'^'^(r) = ^'^•'^P^^ . We remark that if we consider the total of all explicitly-kept 
terms in the integrand of ()3.22|) . the neglected terms would at most be 0{^\ ^^is happens 
because the leading, i.e., order-|^ terms coming from the four pieces denoted 0(pr) in ()3.22|) 
(which are given in terms of 62(^5 ^)) necessarily cancel, as the leading terms of the potential 
V match those of y'^'^'^^ , The above expression is fully equivalent to that found using the 



1 + 



(3.22) 



2nd-order radial WKB approximation for phase shifts in Refs. 



the WKB approximation was used with the Langer-modified radial potential 



[Actually, in Refs 
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but this is inessential for large partial-wave contributions as the difference corresponds to a 
trivial rearrangement of our y-expansion series]. It is not difficult to see that this equivalence 
between the result based on the radial WKB approximation and our present approach using 
the j-expansion persists to even higher orders also. But our new j-expansion for the radial 
proper-time Green function is considerably simpler. 

The /-sum in fl3.22|) can be performed with the help of the Euler-Maclaurin summation 
formula 

°o /-oo 1 1 

Y: m = 2j^ dl f{l) - -f{L) - -f\L) + .... (3.23) 

All terms in this expansion, including the integral term, can be computed analytically. 
The result of this calculation, which is rather lengthy, is given in Appendix O We thus 
obtain an explicit double-integral representation for T^j^j^{A;m) (see the expressions given 
in ()(^9|) - (|(^1H) ). Then the integration over the proper-time variable s can be performed in a 
straightforward manner. After carrying out these s-integrations, we find that the quantity 
rj^j^(A; m) for sufficiently large L is given explicitly by the form 
1 



J>Jl 



{A]m) 



5e Jo 



dr r 



+ / dr 
Jo 



4/i(r)2 + (2/(r) + r/'(r))2 
3 ^ /4L 



Lr 



2r{L + 2)/i(r)V4L2 + mV2 - y {21n (^— ) + 7| {4/i(r)2 + (2/(r) + r/'(r))2} 

{/i(r)(24mV2 + Uh{r)r^ - 75) + 39f{rfr^ + {6r^f\rf + f"{r))r'^ 

1 



O 



L 



Oie), 



(3.24) 



12V4L2 + mV2 

+/(r)(-2/"(r)r^ + 2Af'{ry - 9)} 

where h{r) = /(r)[r^/(r) — 1]. 

For the expression of rj>j^(A;m) we must subtract the renormalization counterterm 
from the expression ()3.24|) . Using ()3.19p . the renormalization counterterm appearing in the 
definition of the renormalized effective action ()2.8p is 

In this radial ansatz case, F^,^ is equal to 



(3.25) 



= -2r],,Jir){r'fir) - l)r'^ + ^{x^r],,, - x,v,Xa){2rf{rf + /'(r))r^ (3.26) 



Hence the counterterm reads 



-(^--7-ln/.^ 



dr r 



4rV(^) - 8r^f{rY + 8f{ry + 4r/' (r)/(r) + rV'(r) 



(3.27) 
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Comparing the expression ()3.27p with ()3.24|) . we see that the divergence terms as e — 
match precisely between the two. This also verifies that our renormalization procedure is 
a gauge-invariant one. The full large partial-wave contribution to the effective action, the 
sum of the expression in (j3.!24|) and minus the result in (j3.27j) . now becomes 

,3 



Tj^j,{A;m)= f dr ^ In {4%^ + (2/(r) + r/'W)^} 

-2r(L + 2)/i(r)V4L2 + m2r2 + |/i(r)(24mV + Uh{r)r^ - 75) 

12 V 4-^2 + m^r^ 

+39/(r)V' + (6rV'(r)' + f"ir)y + /(r)(-2/"(r)/ + 24/'(r)r=^ - 9)} 
+0 (y) . (3.28) 



This generalizes the result of Refs. Q] where the calculation was performed assuming 
the special form /(r) = ,^,2^p2 (i-e., the single instanton solution). 

For any given radial function /(r), one can then consider the sum of this analytic expres- 
sion ()3.28|) for the large partial-wave contribution, and the numerically determined result 
for T j^jj^{A; m) (based on ()2.28|) and the relation ()2.32|) ) to determine the corresponding 
full renormalized effective action. Each has quadratic, linear and log divergences for large 
L. For large L, the L-dependence in the two expressions cancels, as was originally found in 
[1^ . In fact, to improve the numerical efficiency of the full effective action calculation, one 
may extend the large L expression for Tj^j^{A;m) in ()3.28|) . to include also terms up to 
0(-^). This can be done with the help of our y-expansion, using the explicit expression for 
b2{r,r;t) in ^J^. 

2. Case 2 

With ?ij2 given in ()2.16|) . we have the radial potential 

^(M3,*3)(0 = ^ J ' + ^9{r)hH + r^g{r)\ (3.29) 

If we represent F{s) by the form ()2.30b|) and use the correspondence ()3.19|) . the first part 
of ()2.29|) can be expressed by 



/•oo /"OO 

r}>jj^;m) = dr ds 



1 2 
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5: (2/ + 1) ^ Y: {A(,,,3,3)(r,r;.)-A;^f(r,r;.)} 
h=-lt3=±^ 



l=L+\ 



(3.30) 



For the function A(; ;3_t3)(r, r; s) (or A|^^''(r, r; s)) on the right hand side, we may use the 
y-expansion result in (j3.16p with V{r) taken to be equal to the radial potential in (|3.29p 
(the radial potential Vl^f^^ir) = '^'(^+^)+4 y fjj-^g /3-sum and ta-sum can be done explicitly, 
using the formulas : 



E E 



(3.31a) 



+ (/ + l)e 



i(l+l)g{r)s 



(/ + l)e-^('-^)^('^)^}. (3.31b) 



We then perform the /-sum with the help of the Euler-Maclaurin summation formula; this 
produces a double-integral representation in which the integration over the proper-time 
variable s can be executed without too much difficulty. The result of these manipulations is 
1 r 

rj^j^{A; H = ^ dr [8g{rf + 4rg'{r)g{r) + r^g\rf 



+ / dr 
Jo 



90V4L2 + mV2 
1 



+0 



L 



g{rf^AL'^ + m?r^ - ^ |2 In (^—^ + 7| [Sg{rf + Arg\r)g{r) + r'^g{rf] 
[Qr^g{r-f + (20 - 2A<dL)g{rf - 5r{rg"{r) - 12g'{r))g{r) + 15rV(0^} 
0(e). (3.32) 



The field strength appropriate to this Case 2 reads 



Ob \ 

= -2ri^^ag{r)T^ H {x^rj^xs - Xur]^\^)g' {r)T^ , 



(3.33) 



and hence the renormalization counterterm ()3.25|) is given by 
1 /I 



, 00 

2 \ / j„ „3 



7 — In u / dr r 
24 Ve ' ^ J Jo 



ig{rY + Arg'{r)g{r) + r'^g'ij'Y 



(3.34) 



Again we see that the divergent terms as e ^ in (j3.32|) match precisely those of the 
renormalization counterterm. The full large-partial wave contribution, including the renor- 
malization term, is thus given by 

.3 



j>j^(A;m) = / dr 




— In , 

12 \ALJ 



jf) {^9{rf + Arg'{r)g{r) + rV(r)'} 



20 



90\/4L2 + mV2 
-5r{rg"{r) — 12g'{r))g{r) + 15r^5f'(r)^| 



{6r^^(r)^ + (20 - 24:0L'^)g{rf 



(3.35) 



This expression can now be combined with the numerical low partial wave contribution to 
determine the finite renormalized effective action in a radial Yang-Mills background of the 
form Af,{x) = 2ri^y^Xug{r)^. 



IV. CASES WITH ASYMPTOTICALLY UNIFORM FIELD STRENGTHS 

In Ref. the basis for the validity of the y-expansion structure ()3.5|) was a perturbative 
argument: i.e., in the effective potential V{r) = li(^±il±i + V(r), V(r) can be treated 
perturbation to the centrifugal potential term. We then saw in the previous section that the 
resulting explicit expression in ()3.16p . if used to calculate the large partial- wave contribution 
of the effective action, yields a result completely equivalent to that obtained from using the 



quantum mechanical WKB approximation with scattering phase shifts |12L Il6l |. It is not 
immediately clear what happens if the potential is unbounded as r — cxd, in which case in 
the WKB language one should use a bound state analysis rather than a scattering analysis. 
This is precisely the potential form that arises when the gauge field strength approaches 
a nonzero constant value at infinity. Indeed, in our Case 2 with g(r) = — -j (leading to 
uniform self-dual field strengths), we have the quadratic potential (see ()3.29|) ) 

V(M3,*3)W = -45/3^3 + ^r2 . (4.1) 

Thus, the eigenstates of the corresponding radial operator Tij^ consist of only bound states. 
In principle, for the large partial-wave contribution to the effective action, one might still 
try to use the quantum mechanical WKB approximation; but, the WKB approximation for 
bound states has a rather different structure from that for scattering states. 

It is thus an important issue to know whether or not the y-expansion for the radial proper- 
time Green function, considered in the previous section, retains its validity even when the 
potential V(r) blows up for large r, as above. In fact, our y-expansion-based formula in 
(I3.1(i|l does describe the correct asymptotic form valid for all r G (0, oo) even with such 
an unbounded potential; with the proviso that in ()3.16|) . the exponential prefactor e"*^'-^^ 
cannot be replaced by its (truncated) power series in s, for the form to be valid even for very 
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large r. That is, we also have in/rared physics captured correctly by our large-/ limit form of 
the radial proper-time Green function. We demonstrate this below through the treatment of 
an important special case, that of uniform self-dual field strengths for which the quadratic 
potential (14.11) is relevant. Note that the effective action in a uniform self-dual field strength 

pnn 

background has been studied by many authors before |7|, |33|, |3J| , but not using the partial- 
wave-based proper-time formalism. This analysis will also provide a useful comparison when 
we consider, in the sequel to this paper, the effective action in the inhomogeneous background 
described by the form ()2.9p with /(r) = and g{r) = — ^ tanh(r/ro). 

In our procedure the calculation of the one-loop effective action in the gauge field back- 
ground Afj_{x.) = —Tj^yaX^^'T"') starts from the study of the radial proper-time Green function 
A(/^/3^j3)(r, r'; s), satisfying the equation 



d, - dl + ^^^^ ^ ^ - 45/3^3 + ) A(,,,3,3)(r,r'; s) =0, > 0). (4.2) 



Just like the free radial proper-time Green function discussed in Appendix |XJ in this case 
it is possible to find the corresponding radial Green function in a closed form, by using 
(for example) the method of quantum canonical transformations 
expression involving the modified Bessel function 



3511361. It is given by an 



A(M3.3)(r,r'; .) = /,"\ e-^^-"^^^-^^-^+-'^^+^^'-^--/2.+i ( , ^["'^ A ■ (4.3) 
2smh(i5s) \2siYm[Bs) J 

We check the large-/ limit of the expression ()4.3p directly against our general formula ()3.16p . 
as the potential is specialized to the form V{r) = ^'^''^2'*^^ — ABl^t^ + ^r^. As we explained 
earlier, we want our large-/ limit expression to have validity for s satisfying the condition 
< s/^ ^ 0{1). But no restriction follows on the range of the radial coordinate r, and so, 
for the given potential (which diverges quadratically for large r), we want our large-/ limit 
form to be faithful with the true behavior for arbitrarily large value of r^s. With this point 
kept in mind, we may set s = t//^ and use the uniform asymptotic expansion for large orders 
for the modified Bessel function in Appendix [X] (see ()A8|) ) with our expression ()4.3p . in the 
coincident limit, to obtain the large-/ limit form 

t\ Br e ^ ' ' \ vT+i2+i, 



r, r; — 



2sinh(f) V2^(l + 52)1/4 

3x — 5x^ ^ / 1 



l + ^TTT— + 0(7:^)), (4.4) 
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where u = 21 + 1, z = - — ^ttbt:, and x = . ^, . 

For I very large and < t ^ 0{1) (but with no restriction on the value of r), it is possible 
to simplify the complicated exponential term in ()4.4|) by keeping in its exponent only the 
piece — ^ — ^B'^r'^p, namely the leading terms of the given exponent when / becomes very 
large (but for an unrestricted value of r), and expand all the other terms in increasing powers 
of y. Since z is 0{l), this also implies expansion in powers of 4. Then, from (j4.4j) . we obtain 
the following expansion 

^(M3,t3) (r^r;^ 

1 - 37 + hor + + ^Bhht -—An+Oi-]). (4.5) 



It should be noted that, in the same limit, the leading terms of our potential, V{r) = 
^'^^^2^"*"^ ~ ^Bl^t^ + ^r^, are given by ^ + -^r^. Now, according to the same kind of 
reasoning as discussed after ()3.1fi|) . we may replace ()4.5j) by another expansion in which 
the exponential factor at front is assumed by e^^*-^-**^; this rearrangement is harmless for 
arbitrarily large values of here. The result of this rearrangement is to make ()4.5|) turn 
into the structure predicted by our y-expansion formula ()3.16|) . with all factors precisely 
equal. 

We will now show that our radial proper-time Green function ()4.3|1 can be used to red- 



erive the known expression for the effective action. In a uniform self-dual field strength 

nn 

background, it is known that j33|, |3J| 



[Bs] 
(47rs)2 sinh^(5s) 



tr(xs|x) = J,,/ (4.6) 



Thus we have, as this form is used with ()2.3|1 . ()2.4|1 and ()2.8j) . 



(47rs)2 [sinh^(5s) ^ ' 3 
Here, in view of ()2.4p and ()2.22|) . it will suffice to show that 



(4.7) 



E E E (2^ + l)^(M3,t3)(^,^;^)=2vrVtr(x.|x) (4.8) 



i=0,i,l,-^3 = -'t3=±i 
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since rf^x = 2TT^r^dr after the angular integration. Inserting the form ()4.3p for A(i ^3 t3)(r, r; s) 
into the left hand side of ()4.8|) and carrying out the and sums yields the expression 

sinh2(5s) ' ' \^2sinh(5s); 

00 00 2A;+t/ 



sinh'(5s) ^1^0 ^ 'T{k + u + l)k\ 

where we set z/ = 2/ + 1 and z = j^^^, aside from using the power series representation 
of the modified Bessel function. If we change the summation over [k, u) to those over 
{k, n = 2k + v) and use the relation 

. w o ^ z^^^" ^ ^l^{n-2k)sin\i{{n-2k)Bs) 

;tt^o T{k + u + l)k\ ^0 T{n-k + l)k\ 

_^ „2"-isinh(5s)cosh'^-i(Ss) 



1 °o I |i?r^ coth(_Bs) 



n-l 



4 

(here fcm = (f ) if is odd (even)), we then see that the expression in ()4.9j) reduces to 



-^TTTIV (4-11) 



4sinr(5s) 

This coincides with the expression for 2Ti'^r^ tr(xs|x) when (j4.6|) is used for tr(xs|x). Hence 
we have the relation ()4.8p established. 



V. CONCLUSIONS 



In this work we have simplified significantly the calculational method for the one-loop 
effective action developed in Refs. so that any radially symmetric background case 

may now be studied with calculational efficiency. The computation is split into two parts: 
the contribution from low partial waves is calculated numerically using the Gel'fand-Yaglom 
technique, and the contribution from high partial waves has been computed analytically 
using a modified DeWitt expansion. It is no longer necessary to invoke the results of higher- 
order quantum-mechanical WKB approximation explicitly — this is now automatically ac- 
counted for by using the y-expansion for the radial proper-time Green function. The main 
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results are contained in the expressions ()3.28p and ()3.35|) for the analytic behavior of the 
large partial wave contribution to the renormalized effective action, for two general classes of 
radially symmetric gauge fields. Our approach observes gauge invariance, and can be used 
for any mass value for the associated quantum fluctuations. It can also be applied to the 
case with nonvanishing asymptotic backgrounds. In the sequel, we shall report an extensive 
analysis of the Yang-Mills one-loop effective action (not only for scalar matter but also for 
fermion fields as well), taking the radial gauge-field background form of the present work. 
We can then use these results to check for instance the range of validity of the derivative 
expansion Q|- 

In this work we have used the j-expansion to calculate the large partial- wave contribution 
to the effective action. This expansion could alternatively be used to calculate approximately 
the lower partial-wave contributions as well. Aside from the Langer modification js^ which 
can easily be incorporated in our j-expansion, this is effectively what we have done in Ref. 
jlfij l with the Yang-Mills instanton background; there, the instanton determinant was found 
to be good to 5% accuracy. One might be somewhat surprised by this success. But it need 
not be so surprising; observe that the y-expansion as given in (jH.lfij) also serves to generate a 
systematic derivative expansion. In fact, using the expansion ()3.16|) . we have studied several 
cases of one-dimensional functional determinants (including power-like potentials and the 
case of V{x) oc sech^x), to find that the deviation from the exact value is typically not more 
than 5%. This observation could potentially be used to obtain simple approximate estimates 
for general radial background fields. 
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APPENDIX A: THE FREE RADIAL PROPER-TIME GREEN FUNCTION IN 
GENERAL SPACETIME DIMENSION 



In this appendix we find the exphcit form of the free radial proper-time Green's function 
in n spacetime dimension and then discuss its large angular-momentum limit (to facilitate 
the application of our approach in problems with spacetime dimension not equal to four). In 
n dimensions, the consideration of the Laplacian d^j_d^ in generalized spherical coordinates 
leads to the radial differential operator 

a% = — + ^ ^ _ K{K + n-2) 



where r = \lxl + ■ ■ ■ + x^, and n = 0, 1, 2, ■ ■ [From this form, our expression ()2.15|) is 
recovered upon setting n = 4 and k, = 21]. Then, noting that d'^x = r"~^(i"~^f2, we require 
the free radial proper-time Green's function A|5j2(r, r'; s) to satisfy the conditions 

(d (9^ n — ld K(K, + n — 2)] , ^^^^ , , , /. x 

1^ - ^ - -^Yr + I " ' = °' ' > °) ^^^^^ 

s-^0+: Af7-(r, r'; s) ^ ^5(r - r'). (A2b) 
We introduce the modified radial proper-time Green's function A^^^i^, r'; s) according to 



A'^%ir,r';s)=r^A'^%ir,r';sy^. (A3) 



In terms of this function, (0^ and (IX2HI1 can be rewritten as 

- 1^ + V^(':nV0} A^/Kr,^'; s) = 0, (for . > 0) (A4a) 
s^0+: Al72^(r, r'; s) ^ 5(r - r'), (A4b) 
where the centrifugal potential is 



= ^ + „ - 2) H- '" -^f -^' 1 i - f^. (A5) 



To obtain the explicit form of AJ.72) can resort to a variety of methods (developed to 
find the Green function of the one-dimensional Schrodinger equation especially). A particu- 
larly eleg ant method is the one utilizing quantum canonical transformations, as detailed in 
Ref. HQ- 

As it turns out, for A^72) have a simple closed-form expression involving 
the modified Bessel function: 



A^71(r, r'- s) = ^e-i(^^-'^)4,._, i^-j . (A6) 
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Since Iu{z) ~ vfet^ ^^Izf)] large \z\, the s 0+ limit of this expression i 



IS 



s^0+: Al72^(r,r';.)— .-^e-i('--')'{l+0(s)}. (AT) 

The large-K hmiting form of A^''j2(r, r'; s) is of interest. Then, due to the large centrifugal 
potential term in ()A5|) . we expect that the function A^''j2(?", t'; s) be significant (i.e., acquire 
not-too-small amplitude) only when s lies in the range < s ^ A denoting a constant 
of 0(1). Now, for some large given value of n, suppose that we wish to obtain a systematic 
approximation of A^''j2 (?",?"'; s) which can be used for any s satisfying the condition < 
^ 0{1) (this incidentally implies that sn <^ 1). Then, to study the expression in (IA6|1 . 
we use the known large-order asymptotic expansion of the modified Bessel function |37i] 



V large : lv{z) ~ 



V& {l + (z/l/)2}l/4 




(A8) 



with z = and = k + | — 1. [Note that the expansion ()A8jl holds uniformly with respect 
to z (i.e., for any small or large z), and in the limit \z\ oo (for fixed u) goes back to the 
asymptotic form given earlier]. Since we are interested in the case < sk,"^ ^ 0{1), we may 



further take the limit |^| = 1 2s{k+--'i-) I ~^ °° with the the formula ()A8|) (i.e., consider an 
expansion in powers of |^|) and then use it in ()A6|) . After some straightforward calculations, 
we then obtain the large-K expansion of the form 

/ I 4rr' 

We remark that the form ()A9|1 may be used to evaluate a certain quantity which involves, 
say, the integration of A^^^J^i''^, r'; s) over the full s-range (i.e., over s G (0, oo)), as long as k, 
is constrained to be large. This is because, when k is large, (i) A^^j2{''^,r'; s) becomes very 
small unless sk"^ ^ 0(1) (this is also manifest in our form ()A9|1 ) and (ii) for s satisfying the 
condition < sk"^ ^ 0(1) we can exploit the expansion of the form ()A9j) for A^^j2{^,r'] s). 
In view of this, the same purpose can be served by rewriting ()A9|1 as 
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where g'^ = k{k + n - 2) + ("-iK"-3) = + | _ i)2 _ i j^ggg I^K^). Note that this series is 
fully consistent with our formula ()3.16p if we set V{r) to be equal to ^. 

APPENDIX B: THE i-EXPANSION WITH A MATRIX VALUED POTENTIAL 

The y-expansion of A(r, r;s) given in ()3.16p is valid when the potential V{r) is not 
a matrix-type. In this Appendix we shall find a more general form which can be used 
when V{r) and hence also the Green function A(r, r'; s) are matrix-valued. The coefficient 
matrices bk{r, r'; t) in the j-expansion will now have to satisfy the matrix equations in ()3.7a|) 
and (j3.7bjl . Choosing a new independent variable u (instead of r) by setting r = r' + tu and 
writing bk{r' + tu, r'; t) = bk{u, r'; t), we may recast these equations as 



0(r) : dMu, r'; t) + U{r' + tu)ho{u, r'; t) = 0, (Bla) 

1 

(fc = l,2,3, 



0(1^-^") : dMu,r'-,t) + U{r' + tu)h{%r'-,t) - -dlh-iiu,r'-,t) = 0, (Bib) 



The solution of ()Bla|) . satisfying the boundary condition bo{u,r';t = 0) = 1, is 



bo{u,r';t) = P 



-J^U{r'+tiu)dh 



(B2) 



where -P[- ■ ■] denotes the t-ordering. Setting m = in ()B2|) then gives 

6o(r,r;t)=e~*^M, (B3) 

i.e, our formula ()3.1|) for the coincidence limit of bo{r,r';t) holds even when U{r) is matrix 
valued. 

To solve ()Blb|) for higher-order coefficient bk{u,r';t) {k = 1,2, ■ ■ ■), one may follow the 
steps similar to (j3.8p and ()3.9b|) — rewrite the equations as those for the matrix functions 
bk{u,r';t) which are obtained through multiplying bk{u,r';t) from the left by the inverse of 
the t-ordered exponential matrix in ()B2|) . But what we need here is only the coincidence 
limits, i.e., bk{r, r; t) = bk{u = 0, r' = r; t) for small fc, and for the latter it is actually simpler 
to obtain the desired expressions by considering the m = limits of our differential equations 
flBla|) and ()Blb|) and of their derivative relations js^. Specifically, for bi{u = 0,r,t), we 
have the equation (by setting u = with ()Blb|l ) 

9A(0,r;t) + t/(r)6i(0,r;t)-l \dlbo{u,r-t)] =0, (B4) 
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which can readily be integrated (with the 'initial' condition bi{u,r;t = 0) = 0) only if the 
expression for d?ibo(u,r]t) is known. Then, setting m = in the relations obtained 

L J u=0 

after single and twice differentiations of ()Bla|) ) with respect to u, we have 



dt 


dubo{u, 


r-t) 


+ U{r) 


dubo{u, 


r-t)] 


dt 


dlboiu, 


r-t) 


+ Uir) 


'dlboiu, 


r-t)] 



u=0 ' ' 

M=0 

+2tU'ir) IdMu, r- 1)] + t^U"ir)e-'^^'''> = 0, 

L J u=0 

where the result in ()B3|) has been used. From ()B5|1 it follows that 

'duboiu, r; t)] = e-'^^'^ t dti tie''^^''^U'ir)e'''^^'-\ 

L J u=0 Jo 

Using this result, we can go on to integrate ()B6|) to obtain 



(B5) 
(B6) 

(B7) 



dlboiu, r-,t) 



u=0 



Jo 





U{r) 



(B8) 



Now, by using this result in ()B4|) and integrating the resulting equation, we find the expres- 
sion for the coincidence limit 6i(r, r; t)(= &i(0, r; t)): 

h^ir,r-,t) = e-*^M t dt^^h T rfts t2e*^^(^')^7'(r)e-*^^('^) T rfta t3e*^''^'"^f/'(r)e-*^^M 







dt2 t2e*2^(^)f/"(r)e-*^^('-) 



(B9) 



The desired j-expansion, which generalizes ()3.16|) to the case of a matrix- valued potential, 
follows upon using the results ()B3|) and ()B9|) with ()3.5|) . It has the following structure: 

1 



A(r, r; s) 



-sV(r) 



res 



SI 



X ds^s^e'-'^^''Wir)e"-'^^''^ - ds2 sle'^^^'W" ir)e-''^^'^ ^ + O 



.(BIO) 



APPENDIX C: USE OF THE EULER-MACLAURIN SUMMATION FORMULA 

In this appendix we explain how the Euler-Maclaurin summation formula 

HQ can be 

used to sum the various partial-wave contributions to the effective action. First, we present 
the related mathematical theory. Let fix) be a function with continuous derivatives up to 
order 2m + 2 for x G [a, b], where a and b are integers. Then, for the sum 



fin) ^ fia) + /(a + 1) + ■ ■ • + /(6 - 1) + fib), 

n=a 
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(CI) 



we have the Euler-Maclaurin summation formula 

E f{n) = / f{x)dx + -if (a) + m] + Y: 77^,[f''"'-'\b) - J<''-'\a)] + R^, (C2) 

n=a ^ k=l y'^'^l- 

where Bj are the Bernoulh numbers {B2 = |, -B4 = — 3^, -Bg = ^, ■ ■ ■) and the remainder 



term is 



{b - a)B,m+2 fi2m+2)^0^^ ^^^^ q ^ 



(C3) 



(2m + 2)! 

This formula will be particularly useful to evaluate the sum of slowly varying terms with 
decreasing derivatives. (A nice treatment on this formula, including the derivation, is given 



in Ref. 



As for the partial wave sums in our work, we may use ()C1|) in the form (here / 



2' 



being integers) 

00 

E /(O 



E/ 

n=2L 



n 



f{L) 



u:/(i)-4/<-)-n[;|/(i) 

2 |^°°/(/)rf/-l/(L)-l/(L) + ---, 



+ 



f{L) 



x=2L 



(C4) 



assuming /(oo) = /'(oo) = 0, etc. To deal with various Z-sums appearing in ()3.22|) . we may 
here take 

f(l) = e-<^^^'+^'^+^°\bo + hi + b2l^ + ■■■), (a2>0). (C5) 

Then, to perform the /-integral f{l)dl, we change the integration variable from / to t by 
setting 



^/2a2S 2a2 

This will put the function ()C5|) in a simpler form, i.e.. 



/(/) ^ m = e-'\bo + ht + he + ■■■) 

and the resulting integrals can be done by using the formula 



(C6) 



(C7) 



\ {n-\ 
2 



_T2 



T 



2k 



2 

n — 1 



k=0 



{n = odd integer) 
--1 



erfc(T) + Te-^' E 0^ + 



2k 



fc=0 

(n = even integer 



(C8) 
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where (a)„ = a{a + l)...(a + n — 1) is the Pochhammer symboL Now the /-sums in fl3.22|) 
can be performed exphcitly. If we discard contributions vanishing for sufficiently large L, 
these sums for the first two terms in the right hand side of (|3.22j) read 

oo 

^ (2/ + l)(2/ + 2)A(,,^.=,+ i)(r,r;s) 



^s(3i,+2/M-2.VM^),rfc ( y^(2L + 1 + r'fir))] + —(-1 - 8rV(r) 

r / 32s 



sL^ 8rsL2 2r^sL^ 24s^L^ Ss^L^ 
■ — /(O —f (r) ^/ (r) + — ^ ^^^''^ 

-/ (^) - + 77^ - ^/(^) - 3r2sL/(r)2 - -—f'{r) 



3r 3r^ y Y2r2 6 

8rsL , UrhL 2r^sL 2r^sL 
— ^/ W ^ — f{r)f [r) —f (r) —f[r)f (r) 

-fir) + 5 f i'') :r^fir) 5 / (^) 



3^4 3 ' 3 ' ' ' 3r^ ' 3r 
/(r)/ (r) + / (r) — + ^ , / (r) + 



^6 3?^"^ 9?^^ / I 

= J-i(r,s), (C9) 



^ (2/ + l)(2/ + 2)A(,+ i,^.=,)(r,r;s) 



e^(3^+2/M-2.VM^)erfc f v^(2L + 2 - r'f{r))\ + -^(-1 - 8 

\ r J y 8s^ 32s 

+8rV(r)^) + ^(/(r)(r - Ar'f{r)) - Sr^/W^ - r^'W^) 
lo 

V^ivrs I -5 V 2s 2s j 

( 43L 2sL3 les^L^X / 4 Sr^ r%„, , 



V 12 3r2 3r4 ; \ 3 8 ' ' 3 ' ' 6 
2sL2 sL2 8rsL2 2r2sL2 32s2l4 8s2l4 

, 16s2L%,, , 32s3Ln , /2sL 5sL „ , . ^ 2 r ^2 ^HL ^, . , 
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WrsL Ur^sL Ar'^sL 2r^sL 
+ / {r) — /(r)/ (r) —f (r) —f{r)f (r) 



3r^ 3 3 3r2 •' ^ ' ' 3,- 
7, f[nf (r) — / (r) — ——f (r) + 



3^6 3^3 •< ^ ' 9^8 yj 

= -F2(r,s). (CIO) 

If we set /(r) =0 in these expression, they represent the Z-sums for the last two terms of 
()3.22|) (i.e., for those containing the factors e~'^^(')^''^ and e ''+3)^ ^). In this way we obtain 
the explicit double-integral representation for rj^j^(y4; m) of the form 

2 

rj>ji(^; m) = - dr ds / [j^i(r, s) + J^2(r, s) - J^i(r, s)|/(,.)=o - ^2ir, s)|/(,.)=o ■ 

(Cll) 

The /-sums in (|3.30|) can be performed in a similar manner. 
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